Entanglement induced by tailored environments 
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We analyze a system consisting of two spatially separated quantum objects, here modeled as two 
pseudo-spins, coupled with a mesoscopic environment modeled as a bosonic bath. We show that 
by engineering either the dispersion of the spin boson coupling or the environment dimensionality 
- or both - one can in principle tailor the spatial dependence of the induced entanglement on the 
spatial separation between the two spins. In particular we consider one, two and three dimensional 
reservoirs and we find that while for a two or three dimensional reservoir the induced entanglement 
shows an inverse power law dependence on the spin separation, the induced entanglement becomes 
separation independent for a one dimensional reservoir. 
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Introduction 



Among the various schemes that have been recently analyzed to generate entanglement between quantum subsys- 
tems, particular attention has been devoted to mechanisms based on the the interaction with a so called entanglement 
mediator, i.e. a third system which interacts locally with the subsystems one wants to entangle (see e.gi^^). In- 
terestingly enough the role of entanglement mediator can be played by a large quantum reservoir. For example it 
has been shown that two atoms or two quantum dots interacting with the quantized electromagnetic field become 
steadily entangled even at finite temperature*. At the same time it has been recently shown how cold atoms can be 
used to engineer, with a high degree of flexibility, systems described by spin boson Hamiltonians-. This opens the 
possibility to implement new schemes in which the desired entanglement between two microscopic quantum objects 
can be controlled by a suitable design of a mesoscopic reservoir. In the present paper we will briefly analyze how, by 
engineering either the dispersion of the spin coupling or the dimensionality of the reservoir (or both) one can control 
the spatial dependence of the entanglement induced by a bosonic bath on two spatially separated quantum objects. 



The model 



To illustrate the idea consider a system consisting of two two-level systems, labeled a,/3 located at positions r Q ,r^ 
respectively and interacting with a bosonic scalar quantum field. We describe such a system by the following spin 
boson Hamiltonian 



H = Hq + Hi 



(1) 



with 



Ho = f (02 + 05)+ 2>°!« 

k 

H * = EE [(Mr*K + Ki*i)4) (4 + ,7* 



(2) 
(3) 
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where ojq is the energy difference between ground \g) and excited state |e) (we assume h = 1), a z = \e)(e\ — \g)(g\ , 
(j _|_ = \e){g\ and er_ = \g)(e\ are pseudospin operators and a k and at denote bosonic annihilation and creation operators 
respectively, of environment quanta of energy Wfc. The Afc(rj) are position and frequency dependent coupling constants 
between the atom j and field mode k and contain all the information on the spatial dependence of the quantized field 
modes. We will leave their specific functional form unspecified as long as possible in order to simplify the calculations 
and to obtain expressions which arc as general as possible. Note incidentally that the interaction Hamiltonian Hi 
differs from the one examined inl 

Straightforward time independent perturbation theory^, leads to the following expressions for the the perturbed 
(dressed) ground state of the overall system system and for the interaction energy A 
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I*) = c gg \g a9p )\Q k ) + - Q (H, - A) |*) (4) 

A = (g a gp\H^) (5) 

where \g a 9p}\0k) is the unperturbed ground state with energy Eq — —loq and Q =1 — \g a gf3){9agi3\ is a projector 
acting on the atomic degrees of freedom. The state equation (j4|) is in closed form and can be expanded at the desired 
order of approximation. The coefficient c gg is chosen in order to normalize 

The normalized, "dressed" ground state of the system can be written at second order of approximation in the 
following compact form 

l* (2) ) = c gg\9a9p)\0k) + c ee \e a ep}\O k ) + ^ c e<j,fc|e Q g/3)|lfc) + c.ge,k\9c,e.fi)\lk) 

k k 

+ 2 E c 99: kk ' \9a9/3)\HH') + ^ c ee,fcfe' leo-e^) |l fc lfc/) (6) 



The ground state entanglement 



To quantify the entanglement induced by the reservoir between the two pseudospins we will use the two-tangle, i.e. 
the square of the concurrence. We remind that given the density operator p a p of a bipartite system of two qubits, 
the 2- tangle r a \p is defined as 

T alP (p) = [maxtofr-fc-fo-fc}] 2 , (?) 

where (i = f, ..,4) are the square roots of the eigenvalues (in non-increasing order) of the non-Hermitian operator 
p = p(oy®<jy)p*(o y ®<Ty), <7 y is the y-Pauli operator and p* is the complex conjugate of p, in the eigenbasis of a z ®a z 
operator. For states of the form like ©, the 2-tangle takes the following simple forir& 



r a \p = 4 |(e Qe/3 fc |*( 2 )| 2 = 4 |(e Qe/3 fc | 



Q 



-Hi 



Q 



Eq — H o Eq — Hq 



Hi\g a gfjQk) 



E 



A fe (r Q )A^(r /3 ) 



UJ + LOk 



(8) 



In order to extract information on the dependence of the induced entanglement on the spatial separation between 
the two spins we must now specify the functional form of the coupling constants Afc(rj). As far as the mode structure 
is concerned, we assume the usual plane-wave expansion 



Afcfc) -> A k (ri) = e k - 



V 



(9) 



The coupling factor e k contains spin parameters which we do not need to specify for our purposes. We will simply 
assume that it depends on the field mode frequency spectrum via a power law, i.e. 

e k oc lu£ (10) 
where k is the momentum of the oscillators. When V — > oo the sum over momenta k can be replaced by an integral 



E 



A fc (r a )A*(r /3 ) 
w + W/c 



1 



(2tt) s 
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A ,.B ik r 
e k fc k e 

LO + W k 



dk k 



d -i k 2 »f(kr) 
k + k 



(11) 



where r = — r^, d is the spatial dimension and a linear dispersion relation between lo^ and k = |k| has been 
assumed. The function f(kr) comes from the angular integration J dOk and depends both on the spatial dimension 
d and on the dispersion law of the coupling. In the atom-radiation interaction the coupling has vectorial character 
and an r -6 scaling of the tangle is found 9 . Here for the sake of simplicity we assume a scalar coupling. We examine 
now in some detail the dependence of energy and tangle on the spatial dimension d and on the value of v. 
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d = 3 

In three dimensions the angular integration in gives 

f(kr) = J <m k e lU r = Jdcj) Jd9 sin e iferc ™ 9 = 47r^^, (12) 
Assuming = 1/2, as in the case of the coupling with the electromagnetic field, we obtain 

v^MM a if & ^ d3) 

cjo + r J k + k 

The integral present in (fT3"]) formally diverges. The introduction of a cutoff momentum k c (always present in real 
physical systems) gives 

kc k 2 sin kr 1 

dfc— — = -5- {(kq — fe c ) r cosfe c r + smA; c r+ (14) 

fc + fc r 

+ fcor [— 1 + kor sin fcor (Ci(fcof) — Ci(fcor + fc c r)) +fcor cos fcor (— Si(fcor) + Si(fcor + fee?*))]} 

Neglecting the unphysical fast oscillating terms cosfc c r and sinfc c r and and keeping the leading terms in k r <C 1, the 
tangle scales as r -4 with the distance. Clearly, different couplings lead to different scaling laws. For instance, for flat, 
i.e. frequency independent coupling , v = and the tangle scales as r -2 . 

d= 2 

In two dimensions we have 

f(kr) = J di2 k e ik - r = Jd<j)e lkrcos ' 1 ' = 27rJ (fcr), (15) 
where Jo denotes the Bessel function of order 0. Thus 

v«M K [ dkk ^ J jiM (16) 

i£ — ' w + ^7 7 k + k 

Numerical calculations show that this quantity scales approximately as for v = and t/ = 1/2, and as r~ 3 for 
v=\. 

d= 1 

For a one-dimensional environment we get 

v«M«r«^ (17) 

^ w + J_ x k Q + k 

Assuming 1v — n, where n is an integer, we get 

Z-^ r,i„4-/, 1, / J-„ -I- £■ /„ kit ^ ' . 



w + 7-oo ^0 + k J k Q + k J k- ko 



1JL 



>o ifcr /"CO —ikr 

dk k n - - / dk k n - 



u 



ko + k Jo k — ko 

i.e. the tangle turns out to be independent of the spin separation r. 



(-l) n ink'S e~ lkor (18) 
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conclusions 



In summary we have shown that the interaction of two spatially separated pseudospins with a common environment 
leads to an entanglement whose scaling on the spin spatial separation depends on the dimensionality of the environment 
and on how the coupling constants scale with respect to the frequency of the reservoir modes. This result holds for 
the class of physical systems described by a spin-boson Hamiltonian of the form ([3]) . In particular, for a three and 
two dimensional environment we have found a dependence which scales as r~ n , with different values of n depending 
on the dispersion of the coupling constants. On the contrary, in one dimension the tangle does not depend on the 
distance between the two spins. 

Note that our analysis has been carried on by a perturbative expansion of the ground "dressed" state of the overall 
atoms field system. In other words we have studied the equilibrium state of the overall system. Although a full 
time dependent analysis of the atomic entanglement during the process of approach to equilibrium is out of the 
scope of the present articles we would like to point out the difference between our approach and some of the existing 
literature. Some aspects of the time-dependent entanglement which builds up in the collective decay of two subsystems 
interacting with a common bath has been analyzed 

inio^. 

In particular they both study the time dependent as well 
as the asymptotic entanglement relation with the collective damping rate. In both studies however, consistently with 
the Wigner Veisskopf approach to the decay process, the ground state is not dressed. Here instead we have focussed 
our attention on the role of the bath as an entanglement mediator in the ground state. Of course, in the presence 
of collective interactions with a reservoir, there may exist subradiant states which are effectively decoupled from the 
environment. The entanglement which would characterize such state is out of the scope of this manuscript and cannot 
be analyzed with the techniques used above. 

Acknowledgments 

We acknowledge support from the PRIN 2006 "Quantum noise in mesoscopic systems" and from EUROTECH 
S.p.A. 



1 F. Plastina and G. Falci, Phys. Rev. B 67, 224514 (2003). 

2 A. Blais, A. van den Brink, and A. M. Zagoskin Phys. Rev. Lett. 90, 127901 (2003) 

3 J. Majer, J. M. Chow, J. M. Gambetta, Jens Koch, B. R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, 
A. Wallraff, A. Blais, M. H. Devoret, S. M. Girvin and R. J. Schoelkopf Nature 449, 443, (2007) 

4 D. Braun, Phys. Rev. Lett. 89, 277901 (20 02); D. Braun, Phys. Rev. A 72, 062324 (2005). 

5 D. Jaksch and P. Zoller, Annals of Physics 315, 52, (2005) 

6 J. J. Sakurai, Modern Quantum Mechanics, Addison- Wesley (Reading, Massachusetts) (1994). 

7 W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

8 M. A. Cirone, G. Compagno, G. M. Palma, R. Passante, F. Persico, Europhys. Lett. 78, 30003 (2007). 

9 D. P. Craig, T. Thirunamachandran, Molecular Quantum Electrodynamics, Dover (New York) (1998); 

1(1 F. Benatti, R. Floreanini, and M. Piani, Phys. Rev. Lett. 91, 070402 (2003), F. Benatti and R. Floreanini, Phys. Rev. A 

70, 012112 (2004) 
11 Z. Ficek and R. Tanas Physics Reports 372, 369, (2002) 



